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r 1 ^ ■ Abstract 

^. 

The Davey Stewartson hierarchy will be developed based on a set of three matrix dif- 
ferential operators. These equations will act as evolution equations for different types of 
surface deformation in Euclidean four space. The Weierstrass representation for surfaces will 
be developed and its uniqueness up to gauge transformations will be reviewed. Applications 
of the hierarchy will be given with regard to generating deformations of surfaces, and it will 
be shown that the Willmore functional is preserved under this kind of deformation. 



1. INTRODUCTION. 



Surfaces and the dynamics of surfaces play a very essential role in many areas of classical as well 
as quantum physics. Moreover, as far as the area of classical differential geometry is concerned, 
the theory of the immersion and deformations of surfaces has been the subject of intense research 



Domains of study, such as surface waves, deformation of membranes, dynamics of vortex sheets 
as well as certain problems in the area of hydrodynamics are related to the motion of boundaries 
which separate different regions. In particular, in the area of string theory, the action is related to 
the Polyakov integral over surfaces. Certain special classes of surfaces give important contributions 
to various types of physical quantities appearing in these theories and are of interest to consider. 

Recently, Konopelchenko generalized the Weierstrass formulas to the case of generic surfaces in 
M 3 [2-3]. These formulas can be used to study the global properties of surfaces in IR 3 , as well as the 
integrable deformations of such surfaces. This latter aspect is perhaps one of the more important 
reasons for developing these kinds of techniques for inducing surfaces in higher dimensional spaces 
which include Minkowski type spaces as well as Euclidean spaces. These representations can 
then be used to study not only the geometry of surfaces, but the integrable deformations of such 
surfaces as well. 

The generalization of the Weierstrass formulas to generic surfaces in IR 3 which was proposed 
by Konopelchenko consists of the linear system of Dirac equations 



where ip and f are complex-valued functions of z, z G C and p(z, z) is a real-valued function. The 
derivative operators will be abbreviated to d = d/dz and d = d/dz throughout. Using solutions 
to (1.1), three real-valued functions X 1 (z,z), X 2 (z,z) and X 3 (z,z) are defined by the integrals 
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dip = p<£, 



<9y? = -pij}, 



(1.1) 





(1.2) 
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X 3 = - / (ifxpdzf + ip<pdz'), 



r 



where T is an arbitrary curve in C. Regarding the X' l (z, z) as coordinates in IR 3 , (1.1) and (1.2) 
define a conformal immersion of a surface into M 3 . It has been shown that the generalized Weier- 
strass representation of Konopelchenko has several other links to areas of mathematical physics, 
in particular, it can be directly related to the CP 1 nonlinear sigma model [4]. In fact, solutions of 
one system can be transformed into solutions of the other. Moreover both the generalized Weier- 
strass system and the nonlinear sigma model are completely integrable systems [5]. These types 
of representations are useful because in this context, deformations of surfaces can be discussed in 
a straightforward way. Konopelchenko made an important remark in this regard. If a surface is 
represented locally by the generalized Weierstrass representation, then based on the operator 



under the evolution equation induces a local deformation of a surface. The relevant integrable 
evolution equation considered here is the Novikov-Veselov equation, although other integrable 
systems could be considered. Taimanov [6] showed that these formulas for inducing surfaces in 
M 3 describe all surfaces and that the modified Novikov-Vesolov equation deforms tori into tori 
preserving the Willmore functional. These types of application lend an important role to this 
type of representation. The analogous problems for surfaces in M 4 show that this case is very 
different from the three-dimensional case. The main reason is that for tori in M 4 , each equation of 
the Davey-Stewartson (DS) hierarchy describes not one but infinitely many geometrically different 
soliton deformations, and is linked to the fact that the representation is not unique in this case. 

To put it concisely, a surface in IR 3 is constructed from a single vector or spinor function ip 
which is a lift of the Gauss map into nonvanishing spinors, such that it satisfies a Dirac equation, 
and the lift is defined up to a sign. On the other hand, a surface in M 4 is obtained from two 
spinors ip and ip which form again a lift of the Gauss map. However, in this case the lift is defined 
only up to a gauge transformation [7-10] given by , where / is any smooth function. 

It is the purpose in this paper to introduce the system of equations satisfied by these two spinors 




such that the potential p satisfies a particular evolution equation, the deformation in p obtained 
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in R 4 , and the corresponding equations for inducing the corresponding surface. The DS-hierarchy 
will be introduced and a mechanism for deforming these surfaces will be considered. In particular, 
a recipe will be given for obtaining the relevant DS equations pertaining to the first three elements 
of the hierarchy from a matrix system using symbolic manipulation [11]. Applications of these 
evolution equations will be given for deformation of surfaces, and some results related to the 
Weierstrass representation and its gauge invariance will be explored. Finally, it will be shown that 
the Willmore functional is preserved under this kind of deformation. 

2. REPRESENTATION OF SURFACES AND THEIR DEFORMATIONS 

An extension of the generalized Weierstrass representation to four-dimensional Euclidean space 
is based on a pair of spinors ip and ip whose components can be regarded as two independent 
solutions of the generalized Weierstrass system in M 3 , namely (1.1) and (1.2) [12]. Let the spinor 
functions ip = (^1,^2) an d V 9 = (Vi> V2) be defined in a simply connected domain W C C, 
parametrized by the complex variable z, then the components each satisfy the pair of Dirac 
equations 

Vip = 0, Vip = (2.1) 
where T> and T> are matrix operators which are given by 

Before stating how surfaces can be induced from solutions to system (2.1), the following lemma 
will be useful. 

Lemma 1. (a) The components of the spinors ip and cp which satisfy (2.1) also satisfy the 
following conditions 

<9(y^ 2 ) + <9(y^i) = 0, <9(Vw 2 ) - d{(p^ 2 ) = 0. (2.3) 

(b) The one-forms defined by 

Vk = fkdz + f k dz k = 1,2,3,4, (2.4) 
4 



where the coefficients fk are given by 

h = 1(^2 + <pM, h = \(fPii>2 - <pM, 

(2.5) 

h = 1(^1 + <pM, U = 1(^1 - ^1^2)- 

are closed. 

Proof: (a) Expanding the derivatives in (2.3) and substituting (2.1) we obtain that 

d(ip 2 ip 2 ) + <9(</?iV>i) = (d(f 2 )ip2 + ^2(^2) + {9<pMi + Vi(^i) 

The remaining condition is treated the same way. 

(b) Consider the case in which k = 1, the other cases proceed in a similar way. Upon using 
the results from part (a), we obtain that 

drji = Bfi dz Adz + dfi dz A dz 

% — — — % — 

= 2(^2^2) + d(^iV'i)) dz Adz+ -(d(ip^ 2 ) + d(ipi<4)i))dz Adz = 0. 

The next Proposition follows from these. 

Proposition 1. Let the spinor functions ip and (p be defined in a simply connected domain 
W C C and satisfy the Dirac equations (2.1)-(2.2). Then the one-forms rjk in (2.4) define a surface 
in IR 4 by means of the integrals 

X k (z, z) = X k (0) +J Vk , k — 1, 2, 3, 4. (2.6) 

The integral in (2.6) is taken over any path T in W. By Stokes Theorem and Lemma 1, the 
integral in (2.6) does not depend on the choice of path. X 
The induced metric equals 

e 2a dzdz = (\ipil 2 + |^ 2 | 2 )(|<Pi| 2 + IH 2 ) dzdz = mu 2 dzdz, (2.7) 

where U\ = \ipi\ 2 + \ip 2 \ 2 and u 2 = |y?i| 2 + \(fi 2 \ 2 - The mean curvature vector is obtained by 
calculating 

H = ^X z - Z , (2.8) 



and the norm of the mean curvature vector is related to p which appears in matrices (2.2) through 
the expression 

N = y|H|. (2.9) 

For p = p and ip = ±</?, these formulas reduce to the generalized Weierstrass representation for 
surfaces in M 3 . 

In M 3 , a single spinor function is sufficient to obtain a surface. In this case the spinor is a lift 
of the Gauss mapping into nonvanishing spinors, and is required to satisfy a Dirac equation. In 
M 4 , two spinor functions are required to construct a surface, and these functions will form a lift of 
the Gauss map. In fact, not every lift actually satisfies the Dirac equations (2.1). The lifts which 
do are defined only up to gauge transformations. 

One of the reasons for having this type of formalism available to generate surfaces is that 
deformations of surfaces can be obtained and studied in a rigorous way. This constitutes a very 
useful application of these inducing mechanisms. Deformations of a surface are obtained by de- 
forming the potential function p which appears in matrices (2.2) according to some given evolution 
equation. In particular, we will be interested in considering evolution equations which belong to 
the DS hierarchy, as these will appear out of the methodology in due course. To this end, let 
us begin to generate these evolution equations by introducing a corresponding formalism which 
produces them in a precise way. 

To this end, let the operator L be defined as follows 

-U:)+(o :)• 

We consider deformations of the operator L which take the form of a triple of operators L, A and 
B. These operators satisfy 

Lt + [L, A n ] — B n L = 0. (2.11) 
Here t will play the role of the evolution parameter. Now for any nonzero spinor \1/ we have 

[L, d t - A n ]V + B n m = L(d t - A n )* - (d t - A n )L* + B n m 

= Ld t * - LA n * - (dtL)V - Ld t * + A n m> + B n m 



= (-d t L-[L,A n ]+B n L)V. 
Therefore, we conclude that (2.11) implies that 



[L, d t — A n ] + B n L = 0. (2.12) 
Now if L satisfies (2.11), then the solution of the equation 

m = (2.13) 

is evolved according to the equation 

*t = A n m. (2.14) 
Theorem 1. For the case in which n = 1 such that the matrices Ai and Bi in (2.12) are given 

by 

M=( 6 n, Bl J 8 - 8 77 ), (2.i5) 

y p a I \—p — qo — ol 

then (2.12) is exactly equivalent to the Davey-Stewartson II equations 

d t p = dp + dp, d t q = dq + Bq. (2-16) 

Proof: Let ^ = (ipi,ip 2 ) De an arbitrary two-component spinor, then by matrix operations we 
have 

L*=( 1. (2.17) 



— c?Vi + #2 



Then based on these matrices, we calculate 



-p^i + 9^2 \ I d q \ I -pip! + dip 2 
-dipi + qip 2 j \ p d J \ -dtpi + qip2 



d t (-pi>i + dtp 2 ) ~ <9(-Wi + difa) ~ q(-di/>! + #2) 
dti-dipi + qip 2 ) - p(-ptf>i + dip 2 ) - + # 2 ) 



(2.18) 



as well as 



L(d t - Ai)* 



-3 q I \ d t ip 2 - pxpi - Bip 2 



-p(d t ipi - dipi - qip 2 ) + d(d t ip 2 - pipi - dip 2 ) 
-d(d t ipi - - qi[> 2 ) + q(dtiJ 2 - pipi - 8ip 2 ) 



(2.19) 



(2.20) 



with 

Bi l*=( d ~ d -p-A( -p_^+ d ^ 

\-p-q d-d J y -di/ji + #2 

(8 - d)(-pipi + dip 2 ) - (p + q){-di>\ + #2) 

_(p + q)(-p^ + cty 2 ) + (d - 8){-8^ + qip 2 ) 
Substituting (2.18), (2.19) and (2.20) into (2.12) and simplifying, the top element of the resulting 
matrix reduces to 

ipi(dtp -dp- dp) = 0, 
and the lower element of the matrix reduces to 

i>2(-d t q + dq + 8q) = 0. 

These results are exactly system (2.16). & 

These calculations can be best carried out by means of symbolic manipulation. In fact, for the 
cases n = 2 and n — 3, the basic structure of the matrices will be obtained, and then the rest of 
the proof makes use of this [11]. 

It should be noted that the DS I hierarchy is a related system of nonlinear equations which 
are obtained from the DS II hierarchy by replacing the variables z and z by real- valued variables 
x and y. 

Consider the following reduction of the system (2.16) which is specified by taking 

p=-u, q = u. (2.21) 

It is then seen that system (2.16) is compatible under these substitutions and the pair reduces to 
the single expression 

d t u = du + du. (2.22) 
Theorem 2. For the case in which n — 2, the matrices 



— d 2 — v x qd — dq 
—pd + dp d 2 + v 2 



(2.23) 



and 



, <)' 2 ■ <) 2 ■ /•; • r, -(p + q)d + dq - 2dp . 
1 (p + q)d - dp + 2dq -id 2 ■ i) 2 ) - (/'; • r 2 ) i 
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where V\ and v 2 satisfy 

Bv x = -2d(pq), dv 2 = -2d(pq) (2.25) 
generate the following system of equations under (2.12) 

d t p = d 2 p + B 2 p + ( Vl + v 2 )p, d t q = -d 2 q - B 2 q - ( Vl + v 2 )q. (2.26) 

Proof: With L^> given by (2.17), we have that 

(d-A)L*=( 9timi + 92{mi + Mmi ~ qB(m2 + fitiW* 2 ) (2 27) 

{ 1 2 ' \ d t (L^) 2 +pd(L^) 1 -dp(L^) 1 -d 2 (L^) 2 -v 2 (L^) 2 J 

and 

{dt _ = ( W + d 2 ^ + - qdfr + 9# 2 \ = ( A, \ 

\ dti) 2 + pdi)i - dp%pi - d 2 tp 2 - v 2 i) 2 J y A 2 J 

Substituting (2.27), (2.28) and (2.29) into (2.12), the pair of equations in (2.26) is generated by 
calculation, with the p equation the upper component and the q equation as the lower component. 

If we consider the reduction given in (2.21), the q equation in (2.26) is not compatible with 
the p equation. However, by modifying the pair of matrices in a straightforward way, it is possible 
to obtain a compatible pair from (2.12). The following result formalizes this objective. 

Corollary 1. If the matrices (A 2 , B 2 ) in Theorem 2 are replaced by the matrices (iA 2 ,iB 2 ), 
where A 2 and B 2 are given by (2.23) and (2.24) respectively, then (2.12) generates the following 
pair of equations 

p t = i(d 2 p + B 2 p + (vi + v 2 )p), q t = -i{d 2 q + B 2 q + (vi + v 2 )q). (2.30) 

Given Corollary 1, equations (2.30) are compatible under the reduction given in (2.21), which 
reduces the pair of equations given in (2.30) to the single equation 

d t u = i(d 2 u + B 2 u + (v + v)u), Bv = d(\u\ 2 ). (2.31) 

The last system in the hierarchy to be considered here is presented in the next result. 



Theorem 3. For the case in which n — 3, define the matrices 

/ <9 3 + \v 2 d — 3u»i g<9 2 — (Bq)B + <9 2 g + |v 2 g 
I pd 2 — (dp)d + d 2 p + <9 3 + §-u 2 <9 — 3u> 2 

and 




(2.33) 



where 

3 

&n = -^22 = <9 3 - <9 3 - -(t>i<9 - v 2 <9) + 3(w 1 - w 2 ), 
612 = -(p + <?)<9 2 - |(P + <?)^2 - (3<9p - Bq)B - (3d 2 p + B 2 q), (2.34) 

b 21 = -(p + q)8 2 - -{p + - {3dq - dp)d - {3d 2 q + d 2 p). 
z 

such that the Vi and tUj satisfy 

dvi = —2d(pq), dv 2 = —2d(pq), Bw\ = d{pdq), dw 2 = B(qBp). (2.35) 

Then (2.12) reduces to the following pair of equations in terms of p and q, 

d t p = d 3 p + B 3 p + l(v!dp + v 2 Bp) + 3{w l - w 2 + \dvi)p, 

(2.36) 

d t q = d 3 q + B 3 q + \{vidq + v 2 Bq) - 3{w 1 — w 2 — \dv 2 )q. 
This system can be put in the following form by redefining Vi 

d t p = d 3 p + Bp 3 + 3(v 1 dp + v 2 Bp) — 3(d~ 1 [(qBp) z ] + B^ l [(qdp) z })p ) 

(2.37) 

d t q = d 3 q + B 3 q + 3(v 1 dq + t> 2 <9g) — 3(<9 _1 [(p<9<7)z] + B~ 1 [(pdq) z ])q 

The equations in (2.37) are compatible under (2.21) and cause (2.37) to reduce to the single 
equation 

d t u = d 3 u + B 3 u + 3(vdu + vBu) + 3(w + -u/)-u, 

<9t> = <9(|w| 2 ), Bw = d(udu), dw' = B(uBu). (2.38) 

Frequently, (2.31) and (2.38) are referred to as the DS 2 and DS 3 equations, respectively. In 
fact, (2.38) is also compatible with the additional constraint u — u, and reduces to the modified 
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Novikov-Veselov equation. In fact, since dv = d(u 2 ) and dw = \d(du 2 ) = \d(dv) = \d(dv) it 
follows that w = \dv and w' = \dv. Then (2.38) is given by 

d t u = d 3 u + d 3 u + 3(vdu + vdu) + -(dv + dv)u, dv = d(u 2 ). (2.39) 

Zi 

This is exactly the Novikov-Veselov equation. 

3. GAUSS MAP, SURFACES AND DEFORMATIONS OF SURFACES IN M 4 . 

The Grassmannian of oriented two-planes in IR n is modeled by the quadric Q n ~2 in CP n ~ l 
defined by the equation 

n 

E z l = °- 

k=l 

For any point P = (zi, • • • , z n ) on Q n -2 if z k = a k + ib k , we obtain a pair of real vectors A = 
(di, • • • , a n ) and B = (6i, • • • , b n ) which satisfy 

\A\ = \B\, A - B = 0. (3.2) 

These equations are equivalent to (3.1), and A and B cannot be zero since homogeneous coordi- 
nates (zi, ■ ■ ■ ,z n ) of a point in CP 71 ' 1 are not all zero. This pair A, B form an orthogonal basis 
of an oriented two plane II. In M 4 , we give some facts concerning Q2 C CP 3 . There are some 
useful connections between algebraic and differential geometry which are outlined here [13,14]. 
The map o from C x C given by 

a(w 1 , w 2 ) — > (1 + wiw 2 , i(l - wiw 2 ), wi - w 2 , -i(wi + w 2 )), (3.3) 

has the property that a 2 = J2t=i a k = 0- Hence [a] takes values in Q 2 = {[Z] e CP 3 \Z 2 = 0}. 
Moreover, on a(C x C), the related mapping o~ x is given by 

[Z!, z 2 , z 3 , z 4 ) -»• (w!, w 2 ) = ( — , : — )• (3.4) 

Z\ — IZ 2 Z\ — IZ 2 

Thus, [a] is a biholomorphic map from C x C into Q2, which when (u>i, W2) G C x C are considered 
as homogeneous coordinates on CP 1 x CP 1 . It extends to a biholomorphic map of CP 1 x CP 1 
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onto Q 2 . If CP 3 is considered under the Fubini-Study metric of constant holomorphic curvature 
two, the induced metric on Q 2 , expressed in terms of (wi,w 2 ) has the form 

ds > = <W + . (3.5) 

(1 + |-U7l| 2 ) 2 (1 + \W 2 \ 2 ) 2 

This implies that Q 2 is the product of two standard spheres of constant Gauss curvature of two. 

An oriented two-plane in M. 4 is defined by a positively oriented orthonormal basis e\ = 
(ei ; i, • • • , ei^) and e 2 = (e 2> i,-- - 1 e 2 ^) defined up to rotations. There exists a one-to-one cor- 
respondence between components of the and points of the quadric Q C CP 3 defined by (3.1) 

as 

z\ + z\ + z\ + z\ = 0. (3.6) 

The correspondence is given by z\. = + ie2,k for k — 1, • • • ,4. Consider the following change 
of coordinates from ?/j to ^ 

i 1 1 i 

^1=2^1 + ^2), z 2 = -{y 1 -y 2 ), z 3 = -(2/3 + 2/4), ^4= 2(2/3-2/4)- (3.7) 

The quadric in terms of the yi can then be written 

ym = 2/32/4- 

This establishes a correspondence between the space CP 1 x CP 1 and G 4i2 in the form of a product 

G 4 ,2 = CP 1 x CP 1 . (3.8) 

It is this equivalence which allows us to decompose the Gauss map into a pair of maps which can 
be written as 

G = (G i> ,G v ), G:W^G 4 , 2 QeW ^(Xl(Q),X 2 z (Q),X 3 z (Q),Xt(Q)), 
and and G v can be represented in terms of a pair of spinors 

GV = (^2) e CP 1 , G„ = (^1,^2) e CP 1 . (3.10) 

The actual coordinates of the surface can be written in terms of the {jpi , <fi) presented in (2.6). 

12 



Such a decomposition will however not be unique. In fact, the spinors if) and tp will be defined 
only up to a gauge transformation [7-6], that is, a transformation of the form, 

( tp2 ) ( e / '02 ) ' ( <P2 ) ( e - fy 2 ) ' ^ 
Define the functions « a and r a in terms of ip a and <^ a , as follows 

«i = etyi, ri = e _/ v?i 
/t 2 = e / -^2, T2 = e _/ <^ 2 . 

It is clear that the ratios K1/R2 and Ti/f 2 are left invariant by the gauge transformation given by 
(3.11). 

Suppose, for example, a lift is constructed which is based on the initial pair of functions 
(si)S 2 ) = (e td cost], sin 77). Then a pair of functions (^1,^2) can be sought such that the ratio is 
preserved, namely, ^1/^2 = si/s 2 and (^1,^2) satisfy (2.1). From (3.11), we can write 

■01 = e / si, V2 = e / s 2 . (3.13) 

Let us obtain equations satisfied by the functions in the lift as well as /. Requiring that the ifj a 
in (3.13) satisfy Dirac equations (2.1), this will be the case provided that (0,rj) and / satisfy the 
equations 

d{e* sin 77) + pe^ +td cos 77 = 0, 3(e^ +te cos rf) = pe^ sin 77. (3-14) 
Eliminating the function p from the equations in (3.14), we obtain 

df + 089) cos 2 77 = 0. 

Using (3.14) to obtain df, we obtain an expression for p 

p = —e^~^~ l6 (id9 sin 77 cos 77 + drj) . 

Once the ip a have been fully determined, the conditions (2.3) can be used to calculate the com- 
ponents of the second lift G v , since these components must satisfy the remaining equation in 
(2.1). 

In fact, for a certain class of functions /, the entire system (2.1) may be preserved in form, or 
may be said to be gauge invariant. 
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(3.12) 



Proposition 2. If the gauge function / in gauge transformation (3.11) satisfies Bf = 0, 
then system (2.1) is left invariant under (3.11) provided that the potential function p in (2.2) is 
transformed or gauged according to P = pe^~K 

Proof: Suppose the functions ip a and (p a satisfy system (2.1). Differentiating K\ in (3.12) with 
respect to z, we obtain 

A similar result applies to k 2 - Next differentiating T\ we have 

Bti = -(Bf)e~ f (pi + e~ f B(pi = -e~ f (Bf)(pi + pe~ f Lp 2 = Pr 2 . 

A similar result holds for the function r 2 . & 

Consider now deformations of surfaces. Konopelchenko [2] introduced the definition of the 
DS deformations of a surface. Integrable deformations of surfaces generated by the Weierstrass 
formulas will be constructed. As mentioned earlier, this is one of the main applications of the 
generalized Weierstrass representation, as it gives a way to construct integrable deformations of 
immersed surfaces. 

Let us start with surfaces in M 3 . Deformations of the functions ip and cp are considered such 
that there are differential operators A n , B n , C n and D n such that 

i>t n = + B n <p, (p tn = C n ip + D n (p. (3.15) 

For given operators, the compatibility condition of (3.15) with (1.1) gives a nonlinear partial 
differential equation for p. Changing the operators on the right of (3.15) generates an infinite 
hierarchy of integrable equations for p. The deformations of ip and ip described by (3.15) generates 
the corresponding deformations of the corresponding coordinates X l (z,z,t). For example, when 
n = 1, the operators in (3.15) can be written as 

Ai = ad, B 1 = 7p, Ci — aq, Dx—^d, (3.16) 

and system (3.15) turns out to be linear. Equations corresponding to higher values of n are 
nonlinear equations, such as the ones presented in Theorems 2 and 3. 
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Now to generate integrable deformations of surfaces immersed in IR 4 , it is assumed that the 
components of both spinor solutions to (2.1) and (2.2) evolve with parameter t according to (3.15) 
under the same operators A n , B n , C n and D n . The compatibility conditions for (2.1) with (3.15) 
fix the dependence of ip and y? as well as p and q on the parameter t, and consequently define the 
deformations of the corresponding surfaces. Thus, the coordinates X' 1 for % — 1, 2, 3, 4 of the given 
surface which are calculated by means of (2.6) are defined in terms of t as well. Specific cases will 
be governed by different reductions of the hierarchy. This can now be put together and stated in 
a precise way. First, let us note that A n in Theorems 1 to 3 depend on two functional parameters 
p and q and for the case in which p = —u and q = u, define A^ = A n . For the alternate case, 
p = —u and q = u, set A~ = A n . In terms of these new operators, the following Proposition can 
be stated. 

Proposition 3: Let surface E be defined by (2.5) and (2.6) for certain initial spinors ip° 
and ip° which satisfy (2.1), and let p(z,z, t) be a deformation of the potential whose evolution is 
described by the reduced equations (2.16), (2.31) or (2.38). Then the equations (2.5) and (2.6) 
combined with a pair of the following equations 

ipt = Afip, ip t = A^tp, 

il>t = iA£i(>, <p t = -iA^<p, (3.17) 

ip t = A^, ip t = Az<p, 

such that ip t =o = ^ and ip t =o = y?° define deformations of the surface which is governed by the 
reduced evolution equations (2.16), (2.31) or (2.38), respectively. 

Proof: The deformation of p in (2.2) is described by (2.11)-(2.12), so the spinors ip and ip°, 
which satisfy system (2.1), are deformed according to (3.17). Thus, for any t, the resulting spinors 
still satisfy the Dirac equations (2.1). From Proposition 1, however, solutions to this system define 
a new related surface S t in M 4 by means of Weierstrass equations (2.4)-(2.6). This process has 
generated a deformed surface S t such that S coincides with the original surface generated by 
spinors (^ , y?°)- & 

Finally several results related to preservation of surface structure under deformation will be 
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given. 

Proposition 4. Let h be any of the functions ipi<pi, 4 ! i<f2, ip2<Pi, ^2^2 which appear in the 
conservation laws (2.3). Then with respect to DS 2 equation (2.31), each of the functionals 

J(h) = J hdzAdz (3.18) 

is conserved with respect to the evolution parameter, d t J = 0. 

Proof: The claim will be shown for the case in which h = tpifii, the others follow similarly. 
To do this, the evolution equations for ipi and <pi are required. These are given by the second pair 
of matrix equations in (3.17) such that Af are based on the A 2 given in (2.23). These equations 
are given by 

4>it = -*[(-d 2 - v)$! + (ud - du)$ 2 ], <p lt = t[-{d 2 + v)(p x + (ud - du)(p 2 \. (3.19) 

Let 



J = J VWi dz A dz, 
then differentiation proceeds through the integral to give 

d t J = / (V'u^i + V'l^it) dz A dz 



= —i J [(-d 2 ipi - vtp 1 )ip 1 + (udip 2 - (du)tp 2 )ip 1 + ipi(d 2 (pi + v(pi) - ipi(udip 2 - (du)<p 2 )] dz A dz 

= — i J [(udip 2 — (du)'ip 2 )<pi — (ud<p 2 — (du)^)^] dz A dz, 

where (3.19) and integration by parts has been used to simplify this. Integrating by parts once 
more, we obtain 

d t J = — i J [uipidfo + ud(<pi^2) ~ ud((p 2 tjji) — u'tpidtp^ dz A dz 

= — i J \2uLf>idip 2 + u^) 2 d(f>i — 2utjjid(p 2 — u^d'fii] dz A dz. 
Finally, replace both u and u in this by the corresponding derivatives from (2.1) to get 

d t J = — i J [2<9</? 2 <9i/>2 + dipidtpi — 2d r i\) 2 d(p 2 — dipid'tpi] dz A dz 



16 



= — 2i J [d(p 2 dil)2 — dil) 2 dip2\ dz A dz 

= —2i J [— y9 2 <9<9^2 + fi2dd4>2} dz A dz = 0. 
Proposition 5. Let W be the Willmore functional denned as 



W 



= J \u\ 2 dzAdz. (3.20) 



Then the DS 2 deformation (2.31) of tori preserves the Willmore functional defined by (3.20). 
Proof: Differentiating W given in (3.20) with respect to t through the integral , we obtain 

d t W = J {u t u + uu t ) dz A dz. 

Replacing the derivatives with respect to t by the DS 2 equation (2.31), this becomes 

d t W = i J ((d 2 u + B 2 u + 2(v + v)u)u - u(d 2 u + d 2 u + 2(v + v)u)) dz A dz 

= i J (ud 2 u — ud 2 u + ud 2 u — ud 2 u) dz A dz. 
Finally integrating this by parts twice, the required result is obtained 

8 t W = i J {ud 2 u - {d 2 u)u + ud 2 u - {d 2 u)u) dzAdz = 0. 
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